Landau-Zener Probability Reviewed 
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We examine the survival probability for neutrino propagation through matter with variable 
density. We present a new method to calculate the level-crossing probability that differs from 
Landau's method by constant factor, which is relevant in the interpretation of neutrino flux from 
supernova explosion. 
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I. INTRODUCTION 

The study of neutrino masses and mixing is one of the 
most interesting issues in particle physics which has also 
considerable impact on astrophysical and cosmological 
problem. Looking for evidence of mixing neutrino flavors 
during its propagation is one method to detect massive 
neutrinos. If neutrinos propagate through matter, mix- 
ing effects can be enhanced. The electrons in the back- 
ground matter induce the mass to the electron neutrino 
trough a charged current. In non- uniform medium den- 
sity changes on the way of neutrinos therefore the mix- 
ing angle changes during propagation and the eigenstates 
of the Hamiltonian are no more eigenstates of propaga- 
tion. Transitions between mass eigenstates can occur. 
The level crossing probability is known as the Landau- 
Zener probability If density changes slowly enough 
those transition can be neglected so the mass eigenstates 
propagates independently, as it does in the vacuum or 
in a uniform medium. This is called the adiabatic condi- 
tion. The solar neutrino conversion is correctly described 
with the adiabatic condition with accuracy of 10 -7 Q. 
If the density changes rapidly like inside supernovas the 
adiabatic condition is not satisfied, then the probability 
of transition between the mass eigenstates becomes rele- 
vant. 

In this paper we focus our attention in the deduction 
of the level crossing probability expanding the temporal 
evolution operator, we found an general expression for 
this probability and we arrived to the usual one taking 
the first term in the perturbation expansion. In section 
|TT]we briefly review the basic elements for describing neu- 
trino oscillations in a medium, the standard classic prob- 
ability is derived from a geometrical picture. In section 
IIIII We develop a perturbation method to find the tem- 
poral evolution which allow us to find the level crossing 
probability. We found that it differs from Landau-Zener 
probability by a factor 7r 2 /4. 



II. FORMALISM 

In the standard model of neutrinos [3| with #13 ~ a 
neutrino state propagating in the matter is assumed to 
be a linear combination of the flavor states \v e ) and \v a ) 



\v{t)) = v e (t)\v B ) + v a {t)\v a ) 

with \v a ) being a determined linear combination of 
and \v r ). 

The two-neutrino system propagating in matter obeys 
the Schrodinger equation 



d u 
dt 



(1) 



with v — (y e ,v a ) T . Using the Pauli spin matrices <x; in 
the ultra relativistic approximation the Hamiltonian can 
be written as [H 



H = (to + Aoa)er — A (e" 



(2) 



where 2Aoa = V2n e GF, m = (m| + mf)/4E and Ao = 
(m| - m\)/AE > 0. 

In the matter basis, v m = e~ l ^' J2 v, the Hamiltonian is 
diagonalized to 



'He' 



P<T3 — 



(m + A a)(To — Ao3 



where 



and 



A = Aqv 7 ! - 2cos26»a + a 2 , 



cot 2(f) = 



cos 29 - 



sin 29 



which give us two eigenvalues 



E 4 



Vo (a ± V 



1 - 2 cos 29a 



(3) 
(4) 

(5) 
(6) 



associated with the effective masses m± = -J 4EE± . 

III. SEMI-CLASSIC PROBABILITY 

Plotting E± with respect to a we find two hyperbolas 
with the asymptotic behavior 



lu± = to + Ao (a ± (a — cos 29) 



(7) 
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The differences between the curves and the asymptotes 
satisfy 



I E 4 



1 — cos 2(f> sin 



\E- 



■ uj+\ 1 + cos 2(f) cos 2 ■ 



(8) 



so the asymptotes u>± — cos 2 4>E± + sin 2 4>E T represent 
the mean value of the squared mass in the flavor states. 
Furthermore the probability of finding the eigenstate in 



P 2 




Pi 



FIG. 1: Evolution of the probability 



The conversion probability P is known the Landau- Zener 
probability 

IV. QUANTUM PROBABILITY 

Now let us calculate the Landau-Zener probability 
from the Schrodinger equation for the neutrino system. 
Terms proportional to ctq in ([3]) contribute only with an 
overall phase physically meaningless, so we can drop it. 
When neutrinos are produced inside matter the mixing 
angle changes if the density is a function of the position. 
The angle (f> depends on time while the neutrino travel- 
ing in matter. The Schrodinger equation in the matter 
eigenstates now read 



that is 



ie -i<f>*2 —e*°*v m = -Aa 3 u n 
at 



!|"m = (<M - Aa 3 )v m = H e V„ 



(11) 



(12) 



From (fT2")) \(f>\ determine the energy transition between 
the two eigenstates and A give the gap between levels. If 



M 

A 



< 1, 



(13) 



some flavor is given by how close the hyperboles are from 
the asymptotes. We can interpret FIG. Q] in the classi- 
cal way. Let us suppose that N electronic neutrinos are 
produced inside matter, classically there is N\ — NP\ 
neutrinos of mass m+ and = NP2 neutrinos of mass 
77i_. When they go into the vacuum they are going de- 
tected like N e = pOJVi 
type. Then the survival probability is 



P2N2 neutrinos of electronic 



P» e ^v e = ^ = P1P1+P2P2 = i(l+cos2#cos20). (9) 

Actually the neutrinos of mass m + travelling trough mat- 
ter can be converted into neutrinos of mass m_ and vice- 
versa because a quantum tunneling effect. The num- 
ber of conversions must be proportional to the difference 
N\ — N2, so when they travel in the vacuum there will 
be iV° = Ni + P(N2 — Nt) neutrinos of mass m+ and 
N$ = N 2 + P(Ni - N 2 ) neutrinos of mass m_, where 
P is the conversion probability. The number of detected 
electronic neutrinos is iV e = Pf-ATf + P2 and the prob- 
ability for detecting a neutrino electronic is now 

P V .^u. = P?Pl+P2P2-P(Pl-P2)(P? -P 2 °) 

1 / 1 + (1 - 2P) cos 20 cos 20). (10) 



the off-diagonal terms of the effective Hamiltonian H e 
can be neglected and the system of equations for the 
eigenstates decouple. This is the condition of adiabatic- 
ity. For non-adiabatic limit we can not decouple the neu- 
trino system. 

The survival probability for electronic neutrino is 



P, 



(*) = xlU{t)x* 



(14) 



where % Q = (cos a, sin a) T are the components of the 
neutrino electronic in the basis of the Hamiltonian eigen- 
states. U(t) is the temporal evolution operator which 
satisfy the Schrodinger equation 



i-U(t) = H e (t)U(t). 
Because of unitary U (t) can be written as 



U (t) = a (t)<To + ia(t) ■ a, 



(15) 



(16) 



where the at are real and J2i = 1 • The probability 
in function of this parameters is written as@ 
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Pu c ^u e (t) = — (1 + cos 20 cos 20(oq — a\ — a\ + a\) + cos2#sin2</>(aia 3 — aoa 2 ) 
+ sin 20 sin 2(f>(a,Q + a 2 — a 2 — a 2 ) + sin 28 cos 2(f>(aia,3 + aoa 2 

I 



(17) 



The a,i coefficients can be found solving equation (|15p If the Hamiltonian does not commute for different times 

we can do perturbation theory splitting the effective 
Hamiltonian in a no perturbed and perturbation parts, 



which can be written in a differential form 
U(t + dt) = (l - iH e (t)dtjU(t), 



(18) H e = H Q + Hi, in our case 



with the condition U(0) = 1. From (fT5|) it is straightfor- 
ward to find that 



U(t) 



N / N t* 

lim Yl h-iH e (kt/N)dt) = H e - iB ° (t)dt (19) 
N ^ cc k=o ^° 



If [H(ti),H(t 2 )] = for any pair (t u t 2 ) trivially 



e -iHdt _ -i f Q H(t)dt 



Hn 



-A03, H\ — 4>(J 2 . 



(20) Assuming (Ho) > (Hi) we expand flTSJ) as 
I 



(21) 



[/(*) = .Y e -*Ho(t-ctt)«« _ iifl ( t _ L-iH (t)dt _ iHl ^ d ^j L-iH {t+dt)dt _ . Hi q + d ^ dt \ 

-i f* £T dt 

e J o 



i / dtie Jo Hi(ti)e Jtl 
Jo 

t2 m ~* f 1 H " dt tj u \ ~ l IT H ° dt tj u ^ - l ft H ° 
dtie J o Hi(ti)e Jt i Hi(t 2 )e Jt ^ 



dt 2 
o Jo 



This equation can be graphically represented as 



which can be parametrized as 



U(t)= 

l 

with the Feynman rules 
_ -i / H dt 



dtj.... (-iHx(tj)). 



Using (|21|) the time evolution operator can be ex- 
pressed as 



\n=0 J j=l J 



io"3 J Adt 



U{t) = (cosA + i(T2sinAe l/3<T3 )e m<T3 , (23) 

where a is a monotonous function of time and /3 depends 
on the time to reach the vacuum. Comparing with Q16p 
we find for the fully averaged probability (jTTJ) , over the 
time of production and detection, is 



(P»^» B ) = \(l + 



cos 20 cos 20(1 — 2 sin 2 A) ) . (24) 



Comparing with (jTUJ) we can see that the probability 
conversion is given by P — sin 2 A, which is the modulo 
squared of the o 2 coefficient in (|2"2"j) . that is 



xe J o 



(22) 



n=0 J j=l 



(25) 
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This is an exact expression for the Landau-Zener proba- 
bility. At lowest order in (j), the Landau-Zener probability 
is 



and 



Plz = 

From ([5]) we obtain 



dti<fi(ti)e 



i f* 1 Adt 
Jo 



( sin 29 



2(l-2acos26» + a 2 )' 



(26) 



(27) 



Considering that the main contribution is near the reso- 
nance region, a — cos 29 and assuming the neutrinos are 
produced above this region we can extend the limits of 
the integral in (f2"6")) over all a = a — cos 29, 



dt\(j){ti)e 



i f' 1 Adt 
Jo 



10 

where 



sin 29e l1 ^ 
a? + sin 2 29 



1(a) = 2A / V5 



sin 2 29—. 



a 



da (28) 



(29) 



The integral in (|28j) has poles in a — ±i sin 29. This 
integral is calculated to give 



4 



(30) 



To find (|29|) we need to know the functional form of a. 
For example assuming a constant we have 

I(i sin 29) = 2A ° sm22g l n (zsin2g), (31) 



(32) 



where 



_ A 



(33) 



is the adiabatic parameter. Probabilities for other den- 
sity distribution can be found in the literature 0- 

In the usual expression for Landau-Zener probability 
Plz — * 1 when 7^0. It seems that (f3"2"j) is not correct 
because at this limit Plz — * 7r 2 /4 for us. But in this 
situation the perturbation approach (|26p is not valid and 
we need to take the expression (1231) . 



V. CONCLUSIONS 

In this paper we have reviewed the Landau-Zener prob- 
ability starting from standard approach and introducing 
a perturbation method to solve the temporal evolution 
operator. We found that our expression differs from the 
standard one by a multiplicative factor 7r 2 /4 ~ 2.6 which 
at the present experimental resolution is irrelevant, but 
in the interpretation of the neutrino flux from supernova 
explosion [8j could be very important correction because 
of the non adiabatic neutrino propagation. 
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